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Recently we introduced an inflationary setup in which the inflaton fields are matrix valued scalar 
fields with a generic quartic potential, M-flation. In this work we study the landscape of various 
inflationary models arising from M-flation. The landscape of the inflationary potential arises from 
the dynamics of concentric multiple branes in appropriate flux compactifications of string theory. 
After discussing the classical landscape of the theory we study the possibility of transition among 
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various inflationary models appearing at different points on the landscape, mapping the quantum 
landscape of M-flation. As specific examples, we study some two-field inflationary models arising 
from this theory in the landscape. 
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I. INTRODUCTION 



The inflationary paradigm, the idea that the early Universe has undergone a nearly 
exponential expansion phase, has appeared as the leading candidate for explaining the recent 
cosmological observations data [1]. The simplest and still successful model of inflation is a 
massive, free scalar field minimally coupled to the Einstein gravity Nonetheless, motivated 
by various beyond the Standard Model particle physics or supergravity and string theory 
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settings, many models of inflation have been constructed by introducing more non-trivial 
potentials for the scalar field and/or the addition of other scalar fields to the model. 

The option of having a large number of scalar fields, instead of a handful of them, has 
been particularly motivated by string theory inspired models, where after fixing the moduli 
in compactifications they appear as light scalar fields subject to different potentials [2j]. 
A large number of (decoupled) scalars has the advantage that the contribution of each of 
them to the Hubble expansion parameter during inflation adds up, leading to a substantial 
number of e-folds, even if the potential for each field is not flat enough to sustain a successful 
period of inflation On the observational sides, a generic multiple field inflationary system 
predicts entropic and non-Gaussian perturbations which both are under intense observational 
constraints [1]. 

in y, again motivated by string and D-branes settings, we introduced a model of inflation 
which has a large number of fields. In this model the inflaton fields were taken to be matrix 
valued objects and hence the model was dubbed as Matrix Inflation, or M-flation in short. 
The non-Abelian and non-commuting nature of matrices play a crucial role in M-flation 
construction. 

pi 

We showed in [4J that in a special corner of the rich "landscape" of the M-flation, one can 
reduce the theory to a standard chaotic inflationary model with at most a quartic potential. 
The matrix nature, however, now shows up in removing (or easing) the unnaturally small 
values of the couplings for these models. In this work we would like to analyze the landscape 
potential of M-flation in more detail. The landscape arises from the dynamics of coincident 
branes subject to appropriate fluxes in string theory compactification. At different points 
of this landscape we have inflationary models which are classically disconnected and the 
parameters at each point could be adjusted to achieve a period of slow-roll inflation with 
enough number of e-folds and acceptable spectral index. These different vacua, in principle, 
can tunnel to each other via quantum effects. If some fine-tunings of the couplings are 
tolerated, similar to the old inflation idea which should have terminated through first order 
phase transition [5], slow-roll inflation could end through bubble collision. 

The paper is organized as follows. After reviewing the basic setup of M-flation in section 
HT| we study the landscape of M-flation theories in section IIHI In section [IV] we analyze 
the landscape of the potential and hence inflationary models for a given Matrix inflation 
theory. Here we give a counting of various possible models (a counting on the landscape). 
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This is basically a counting of number of N x N reducible representations of SU(2), which 
is much larger than N 2 for large N. We also discuss the tunneling between various vacua. 
As specific examples, we study some two-field inflationary models arising in our setup. In 
section |V] we study numerically the physical predictions of these two-field models for CMB 
perturbations. The conclusions and discussions are provided in section I VI I In Appendix 
A we present the cosmological perturbation theory for two field models in some details. In 
Appendix B, we consider the case where we have a positive energy metastable vacuum and 
find out the set of couplings where the first order phase transition from metastable vacuum 
to the true one is possible via quantum tunneling. As stated above, for a small window of 
parameters, nucleation rate is large enough to allow for a a first order phase transition via 
Coleman-De Luccia phase transition. 



II. M-FLATION, THE SETUP 



As in |4| we consider M-flation the inflationary model in which the inflaton fields are 
considered to be matrices and the inflationary potential is constructed from the matrices 
and their commutators. The action is 

8 = J d 4 x^^it!-i^Tr(^^-mJ^^])j , (n.l) 

where the reduced Planck mass is Mp 2 = 8nG with G being the Newton constant and the 
signature of the metric is (— , +, +, +). Here the index % counts the number of matrices and 
we take it to be i = 1, 2, 3. The kinetic energy of has the standard form and $j are 
minimally coupled to gravity. To be specific, as in [4|, we consider the following potential 

V = Tr (-^, $,] + je jkl [*k, *i]*J + ^i) , (H.2) 



which is quadratic in $j or $j]; (1II.2j) is the most general potential with this property. 
We choose the three parameters A, k and m 2 to be non- negative. 

The above action enjoys a global U(N) symmetry under which <3>j are in the adjoint. One 
may also consider the theory in which this U (N) is gauged. The latter is done by replacing 
the partial derivatives of the scalars $j by covariant derivatives, = <9 M $j + ig[Ap, 

where is the gauge field, and by adding the Yang-Mills term for the gauge fields. For 
our purposes, during (slow-roll) inflation and as far as the classical dynamics of the scalars 



is involved, one may consistently set the gauge fields to zero, in which case the inflationary 
dynamics of the gauged and un-gauged models become identical. These two models, however, 
can have a different spectrum of linear perturbations about the inflationary background. 

In 4| we argued that this system, and in particular when the U (N) symmetry is gauged, 
is strongly motivated from string theory where $j represents the low energy dynamics of 
coincident branes in some specific flux compactifications In this picture, N coincident 
D3-branes extended along our Universe, subject to RR six-form field CVg), and through the 
Myers effect 7i], can blow up into D5-branes, wrapping around a two-dimensional sphere in 
the extra three dimensions which are parameterized by matrix valued scalars $j. 1 This two- 
sphere for finite N is a fuzzy two-sphere and in the large N limit it becomes a commutative 
round sphere. In |4| we studied the case where N x N matrices exhibit a single fuzzy two- 
sphere. It is notable that in the large N limit the sector of M-flation studied in |4j , similar 
js], {9] and 10 1, is closely connected to the models of D5-branes wrapping two cycles 



to 



of an internal (Calabi-Yau) space 11 



I2I ]. In these models the D5-branes are moving in 
a Klebanov-Strassler throat 13[ where the DBI effects become important. In this work, 
however, we would like to focus on the matrix effects. As we show below due to the matrix 
nature of the fields we have the possibility of multi-field inflationary models which does not 
exist in the analysis of 
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As we will show in this work, there is also the possibility that we have a large number 
of concentric fuzzy spheres of various radii, corresponding to a picture where we have more 
than one fuzzy sphere or wrapped D5-brane. Moreover, there is a distinct inflationary 
model associated with each of these multi-fuzzy spheres solutions: For the case of multi 



ds 
i 



fuzzy spheres solution we obtain a multiple field inflationary model, where these scalar fie 
are classically decoupled from each other. These fields, similarly to the single field of 
have the geometric interpretation of radii of the fuzzy spheres. To see this we analyze 
the equations of motion of the theory. Starting with an isotropic and homogenous FRW 
background 

ds 2 = -dt 2 + a(t) 2 dx 2 , (II.3) 



1 We note that, as discussed in [4|], the action (III.ll) with the potential (|II.2[) is the action in the lowest 
order in string scale a' . Considering a' corrections adds term higher order in $j. Therefore the string 
theory picture for our model is valid if these a' corrections are small. 



6 



the equation of motions are 

R2= 3^ {~l Tl (^Si^) + V ^ [$ " (IL4a) 
<t t + 3H^ + A [fy, [$,, $,] ] + i $ fe ] + m 2 $, = , (II.4b) 

^ j 

where H = a/a is the Hubble expansion rate. 

It is straightforward to check that the above equations of motion can be solved with $j 
in the matrix form 

$>i = Y,kJT , i = l,2,3 (II.5) 

a 

where J° are iV Q x jV a matrices satisfying the SU{2) algebra 

[J? , J/] = ic^ J fc Q , Tr(J"jf) = \N a {N 2 a - 1)8^ . (11.6) 

Here iV^ are arbitrary non-negative integers subject to the condition N a = N. That 
is, Jf form reducible N x N representation of SU(2). In the above, the value of a 
specifies the number of irreducible SU(2) blocks in the generic N x N matrix. Therefore, 
the range of a can vary from one, corresponding to a single irreducible representation which 
was discussed in some details in (4], to N, where basically we have the solution $j = 0. 

In fact, one can show that ( \ll.5\i is the most general solution to (III. 4b ). Moreover, sum of 
two solutions of the form (III.5[) is not a solution to fill. 4b ). With the above decomposition 
one finds that the equations for cj) a decouple and one may analyze them separately. It is 
then convenient to rewrite the action in terms of (f) a , rather than the matrices $j. Doing so, 
and after rescaling 

n ~r 1/2 

k= - A N a {N 2 a -l) <j> Q , (II.7) 
we obtain a multi-field canonically normalized action for the scalars <p a with the potential 



a 

where 

8A 2 k 

K ~ N a {Nl - 1) ' Ka ~ y/N a (N*-l) [ } 

As we see the fields <p a are decoupled from each other, which is in fact the result (or ad- 
vantage) of using the SU(2) basis we have introduced. Note that although in the potential 



( 1II.8|) there is no interaction between <p Q , they all couple to gravity and contribute to H. 
Therefore, their dynamics are coupled through gravity. 

From now on we call each reducible solution which minimizes the potential a "vacuum" 
solution. In |4J we analyzed the irreducible vacuum, where we studied the N a = N case in 
which a = 1 and as such this vacuum may also be called the "single block vacuum" or the 
"single giant vacuum" (in a reference to the blown-up D3 branes picture). Similarly, when 
a ranges from one to n we have an n-block vacuum or n-giant vacuum. 



III. THE LANDSCAPE OF M-FLATION THEORIES 

There are four parameters in the original M-flation action, A, k, m 2 and the size of 
matrices N. As discussed above, the classical solutions of the theory are described by 
the set of parameters (A a , K a ,m 2 ) and the set of N a which are subject to J2a=i N a — N. 
Depending on parameters (X a , K, a ,m 2 ) one obtains potentials of different forms. Since the 
potential ( \ll.8\i is the sum of potentials which only depend on a single field <p a , we may 
analyze each of them independently. The latter is basically the shape of the potential about 
the single-block vacuum which is studied first. We then analyze the multi-field case. 



A. Analysis of the potential around the single-block (irreducible) vacuum 

As a starter, we briefly review, summarize and expand on the results of J4], where Matrix 
Inflation theory is effectively reduced to a single scalar field theory corresponding to a = 1 
in Eq. ( )H.5p . For that purpose we assume the N x N matrix part of to be proportional to 
the generators of SU (2) algebra in its N x N irreducible representation. Upon the rescaling 

-1/2 



-N(N 2 - 1) 



(HI 1) 



the action reduces to a canonically normalized scalar field theory coupled to Einstein gravity 
with the potential 

VsM = V - + ^ 2 (IH.2) 



where 



A, 



8A 

N(N 2 - 1) 



Ik 



VN(N 2 - 1) 



(III.3) 
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FIG. 1: Different possible shapes of the single block potential Eq. (|III.2|) . 



Depending on the value of the ratio m 2 A^/K 2 = 2m 2 \/n we may have five different cases. 
Note the remarkable property that this ratio is N- independent. In Fig. [T]we have plotted 
the different shapes of the potential. Since in the string theory setup, <f) a represents the 
radius of each fuzzy sphere, we take (fi a to be positive. 

• Case I) k 2 < 2m 2 \. 

In this case the potential has only a single minimum at <fi = for which the potential 
vanishes. All the k = cases, including A0 4 and m 2 2 chaotic inflation potentials are 
in this class. 
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For the rest of > 2m A and hence one may use the parameterization 

2m 2 \ 

sin 2 9 = — 'f- 6 G [0, tt/2]. (III.4) 



K 2 



Case II) 6 = f . 

ot 

of potential at the inflection point is 



In this case the potential has an inflection point which happens at 0o — j^- The value 



,4 



VS = N(N 2 - 1) . 

s 384A 3 V ; 



As discussed in [4j the inflationary model based on this case has a small red spectral 
index and is on the verge of being ruled out. 

Case III) sm-\l) < < f . 

In this case the potential besides the zero energy minimum at = has also a 



minimum at 



and the energy 



where 



* (! -= ^(1 + cosG) , (III.5) 



A, 



V s ° = iV(iV 2 - 1) F(6) , (III.6) 



F(9) = (1 -3cos6)(l + cose) 3 . 

In this case V$ > 0. If the inflationary dynamics happens around the 0o minimum 
such that the end point of the slow-roll inflationary phase is at <fro, then the situation 
is like old inflation and we are back to the "graceful exit" problem. In Appendix B, 
we have determined the set of parameters for which nucleation rate via Coleman-De 
Luccia [ljj quantum tunneling is substantial to cause a first order phase transition 
at <f) Q . We have also shown that it is not possible to tunnel from the false vacuum 
to the true one by Hawking-Moss 15( phase transition, if one demands the COBE 



normalization for density perturbations in the subsequent slow-roll phase. 



Case IV) 6 = sin 



>9> 



In this case 4>o = Ak^/SX^, and F(Q) = so the energy of the potential at the minimum 



vanishes. In [4J this case was called the "symmetry breaking" potential. This case 
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does not suffer from the secondary old inflation phase and hence is a preferred case 
for building an inflationary scenario. Moreover, the potential in this case has the 
remarkable property that it could be obtained from a cubic superpotential and hence 
the model can be embedded in a supersymmetric theory This latter gives us a better 
control over the running of the parameters and the Coleman- Weinberg corrections. In 
string theory set up, the minimum at 0o corresponds to the configuration when N 
coincident D3-branes blow up into the supersymmetric configuration of a giant fuzzy 
sphere with radius determined by 0o- The other vacuum, = 0, corresponds to the 
configuration of commuting matrices with a shrinking size sphere. 

• Case V) < 6 < sin -1 ®. 

In this case F{&) < and the global minimum is at 0o with an AdS vacuum. The 
m 2 = case, which is the case usually studied in the D-brane setting in string theory 
(see e.g. [3, Q, Il6]), falls in this class. It is possible to have slow-roll inflation around 
the minimum = 0. The minimum at 0n is an AdS type and is not suitable as an 
end point for inflation. As discussed in [14 1 the gravitational effects makes the = 
vacuum stable against tunneling to true vacuum at O . 

In the specific examples of two-field models studied in the following section, we restrict our 
analysis mainly to the cases I) and IV) above. 



B. Analysis of the potential for n-block vacua 

We now consider a general n-block vacuum, corresponding to the n-field inflationary 
potential (1II.8P with fjll.9j) . In this case the equation of motion for the fields <p a (when 
gravity is turned off) are essentially decoupled from each other. Nonetheless, all of them are 
coupled to the background metric and give contributions to the Hubble parameter and hence 
the field a will feel the effects of the other fields through the evolution of the background 
metric. The equations of motion for the n-field theory are 




(III.7a) 



a + 3if0 a + d^K = , (III.7b) 
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where V a (4> a ) is given by Eq. (1II.8j) . 

Noting that the angle B given by Eq. flIII.4j) is independent of N a , the form of the 



potentials V a (in the sense that which of the five cases of previous subsection they fall into) 
is independent of a and hence for a generic n-block vacuum we still have the same five cases 
discussed above. The value of the potential at the minimum and at which Q the minimum 
happens, however, depend on N a . 

In summary, irrespective of which set of {N a } is chosen, the generic shape of the infla- 
tionary potentials falls into one of the five cases discussed above and is completely specified 
by the dimensionless ratio k 2 /(m 2 A). 

IV. LANDSCAPE OF n-BLOCK MATRIX INFLATION 

In the previous section we analyzed different possibilities for the generic behavior of the 
inflationary potential depending on how the original parameters of M-flation is chosen. In 
this section we study the landscape of various inflationary models arising from inflation 
assuming a given set of parameters (A, k, m 2 ) and N. 

A. Counting the number of vacua in the landscape 

As the first piece of information on the rich landscape of the potentials, which all are in 
the generic form of ( 1II.8j) with (1II.9j) . we present a counting of these vacua. The case of single 



block, n=l, was studied in IIII Al and depending on the ratio m 2 X/K 2 we may have either 
of the five distinct cases discussed earlier. Next, we may consider the double block vacuum, 
which leads to a two-field inflationary potential. In this case, the number of possibilities 
is the same as the ways one can partition N into two positive integers. Denoting the size 
of two blocks by Aq and N 2 , Aq + N 2 = N, with Aq > N 2 > 1, one has [|iV] possible 
solutions, where [x] is the integer part of x. (Note that exchanging Aq and N 2 does not lead 
to a different theory.) In other words, there are [-AT] possible two-field models coming out 
of M-flation. 

In a similar way one may count Af n (N), the number of independent n- field inflationary 
models arising from M-flation (which is the same as the number of vacua or the minima of 
the potential). J\f n (N) is the number of ways an integer N can be partitioned into n positive 
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integers iVj such that Ya=i Ni = N, 1 < iVj < iV i+ i. One can show that 



Af n (N) = J2^n-i(N + n-l- nk), M 1 = 1, 



(IV. 1) 



fc=i 



for n < y and 



N n {N)=M N . n {2{N -n)) , 



(IV.2) 



for n > y. From the above recursion relations one can compute Af n (N). For example, for 
n — 2, J\f 2 — [y] and for n = N, J\f N = 1. For large A/", and when n N one can show 



that Af n (N) ~ ^(iV - 1) • • • (TV - n + 1) ~ ^t^™ -1 - The number of points (theories) in the 
landscape of the possible theories, for a given set of parameters (A, n,m 2 ) and iV is 



for large N. 

It is worth noting that around each minimum point in this landscape we have an n-field 
model specified with {N a }, where a — 1, 2, • • • , n. Around each point there are 3N 2 —n fields 
where A^ = Yla=i ^ na ^ are no ^ classically excited. These fields, which are not classically 
coupled to the n inflationary fields, can appear through quantum excitations of the fields, 
leading to 3A^ 2 — n isocurvature fields. In summary, among the cosmic perturbations around 
the specific n-block vacuum, we have one usual curvature (adiabatic) mode, n — 1 entropy 
modes and 3A^ 2 — n isocurvature modes. 

B. Transition between multi-field scenarios 

As discussed, if we start with the initial condition that the matrix valued fields $j and 
their time derivative are given in terms of the generic SU (2) (reducible) representation, 
the dynamics of M-flation is such that they always remain in the same sector. In this sense 
the classical (inflationary) dynamics around the "multi giant vacua" decouple from each 
other and one may build an inflationary model around either of these. If we start with 
a field which is initially in the sector specified by a given set of N a , then N a remains a 
conserved quantity by the classical trajectory of the system. In general various fields in 
the same sector specified by a set of N a can mix with each other. That is, in general 




(IV.3) 



71=1 
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the inflationary trajectory in the space of (fi a is curved. We will discuss the details of the 
inflationary dynamics for the special two-field case in the next section. 

This "classical" decoupling can, however, break due to quantum effects and one can find 
various quantum paths connecting vacua with different sets of N a . To see this let us suppose 
that the ratio m 2 \/n 2 takes a generic value and for the ease of the analysis consider the 
epoch the inflationary dynamics has ended and the fields are settled in the minimum of their 
potential. Let us study the tunneling between vacua corresponding to two SU(2) reducible 
representations specified by sets of and which will respectively be denoted by 

(2) 

and Jj . One such path between the two vacua is 

= <f>{t)jP + (1 - 0(t)) Jf ] , 0(t = 0) = 1, <f>(t = t ) = 0. (IV.4) 

One may then evaluate and minimize the action for the above path and use the WKB 
approximation to compute the transition amplitude by exponentiating the value of the Eu- 
clidean action. (We comment that there are other paths over which the action is smaller 
compared to the one evaluated for the above path. In such cases these paths would dominate 
the tunneling.) 

The transition may happen from a local minimum with higher (or equal) energy to a local 
minimum with lower (or the same) energy. To be specific let us compare the irreducible 
(single-giant) vacuum with a generic n-giant vacuum. Depending on the value of the ratio 
m 2 X/ k 2 the single or n-giant vacua could have a higher energy. It is readily seen that 

vs° - e ficld = Z(N?-if Vs ° ' ( IV - 5 ) 

and hence V$ — V®_ md has the same sign as Vg. Explicitly, for the case III the n-giant 
vacuum is more stable and the single giant vacuum can decay (tunnel) into the multi-field 
vacua, as well as tunneling into its own = vacuum. For the case V, we have an opposite 
situation and the single giant vacuum is the most stable one and the multi-giant vacua can 



tunnel into the single giant one. Such an analysis has been carried out in 16|. The symmetry 



breaking IV, is special in the sense that Vg = and hence all the vacua have 

the same energy. The calculation of the transition amplitude for this case has been carried 
out in some detail in [17]. In the analysis there it has been shown that the path which 



minimizes the action is not of the form (1IV.4I) . In 17|] it is shown that the height of the 



potential between the two vacua grows like N 2 (for large N), unlike the N 3 expected from 



a path like (JTVJJ. 

14 



FIG. 2: Two possible paths for quantum tunneling between two block and single block vacua. 
The top figure shows a paths in which the spherical symmetry of the giants are preserved along 
the path. The figure below shows a path in which the transition happens through nucleation of a 
throat between the two giants, through which the inner giant is eventually dissolved into the outer 
one. Both of these paths are coming from quantum fluctuations of the two initial giants and can 
be understood as specific modes of open strings stretched between the two giants. Analysis of 
shows that the path depicted in the lower figure dominates the transition. 

One may also provide a picture for the above tunneling in terms of transition between a 
multi-giants configuration to a single giant one. Let us consider the specific case of transition 
between two and single giant configurations. As was discussed the above reducible solution 
with two blocks (in the large Ni,N 2 limit) has the interpretation of concentric spherical 
D-branes whose radii, A r 1 1 /2 and N 2 4>2/^ are changing in time. Analysis of [l^ suggests a 
transition between the two is dominated by the path which is mediated through nucleation 
of a throat between the two giants, which eventually dissolves them into a single giant (or 
vice versa), as depicted in Fig. [2j These throats are basically the virtual open string loops 
which are stretched between the giants (branes). 

Moreover, there is also a non-zero tunneling amplitude from the non-trivial "giant" gravi- 
ton vacuum to the vacuum at <fi = 0. For the case of supersymmetric theories some of these 
paths may be forbidden due to supersymmetry. 
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V. TWO-BLOCK MATRIX INFLATION 



We have discussed that dynamics in a sector with a given set of N a , a = 1,2 , • • • ,n 
can decouple from the rest of the theory if we initially start in that sector. Although the 
inflationary potential in this sector contains n decoupled scalars, they are in general coupled 
to each other through gravitational interactions, leading to n-field inflationary models. In 
this section we analyze in some details the simplest of such models, i.e. the two-field case 
for which the field <3>j can be expanded as 

*i = 4>(t)j? ) + %(t)4* ) , (v.i) 

where and jf^ are respectively N\ and iVV dimensional (Ni + No = N) irreducible 
representations of 577(2) algebra. In the string theory setup studied in |4|, the fields <fi and 
X indicate the radius of the giant spheres to which the stacks of coincident branes are blown 
up. Since jf } and jf ] are Hermitian, we conclude that and x are rea l scalar fields. 

After rescaling the (j) a , A and k as in (111.71) and (III. 91) we obtain the effective two field 
potential 

Vnfa x) = + ^X 4 - - ^X 3 + ^ (0 2 + X 2 ) ■ (V.2) 

The effective couplings A^, A x , and k x are given by 

with a similar expressions for A x , n x with iVi replaced by N 2 . 

With the initial condition of solutions given by Eq. ( IV. lj) . our system is basically reduced 
into a two-field inflationary system of <fi and x- The background equations of motion are 



H 



3M 2 P 



\<P 2 + \x 2 + v D (<p, x ) 



+ 3i70 + ^y D = O , x + 3Hx + d x V D = 0. (V.4) 

Since m 2 X^/K^ = m 2 X x /K x , depending on the value of this ratio, the potential Vp can 
be classified into exactly the same five cases studied in subsection (IIII A[) and hence we 
do not repeat them here. Below we study some inflationary backgrounds constructed from 
Vd{4>,x)- F° r some details of the treatment of the cosmological perturbation theory for 
two-field models see Appendix [S] 

16 



For general initial conditions, the inflationary trajectory is curved in the (j)~x plane. The 
general result of such an effect will be a nonzero correlation between the adiabatic and the 
entropy perturbations. The magnitude of the effect, however, depends on the values of the 
couplings and the initial conditions. We consider the following cases, which shows to what 
extent the results are initial condition-dependent. 



A. Quartic Potential 

Suppose k = m = 0. The background inflationary potential 

V = + ^X 4 , (V.5) 

has the form of case I. 

Considering number of e-folds, N e , as the clock Hdt = —dN e , in the slow- roll limit one 
can check that 

(f)~ ^y8N~ e Mp sin a , \ - \f&N~ e M P cos a , (V.6) 
where a is an angle. Furthermore, the trajectory in (0, x) space or (N e , a) space is 

o ]\/[2 

^ = R 2 X ~2 + C ; ]V e = ^( s in- 2 a-i? 2 cos- 2 a) (V.7) 
where C is a constant of integration and R is the ratio of the couplings 

2 _ \„ N 2 (iq - 1) 

R = A7 " N l{ Nl - 1) ■ (V ' 8) 



The analysis is similar to 18|, [19( and 20] . 



For the numerical study, we first consider the case where the quartic coupling of x an d ip 
are of the same order. In this case the trajectory is slightly curved, if the initial conditions 
of none of the two fields are set to zero. For definiteness, we assume that 

A = 2.05 x 10~ 13 , A x = 5 x 10~ 13 

<f>i = 22.437 M P , Xi = 8M P , (V.9) 

where the initial value of fields, <pi and x%i are given 71 e-folds before the end of inflation. 
Choosing the natural bare value of A ~ 1, from Eq. (1V.3I) . the above values for the effective 
quartic couplings correspond to the following values for the dimensions of the blocks: 

Ni ~ 34000 , N 2 ~ 25000 . (V.10) 
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FIG. 3: Left graph shows the trajectory for two quartic potentials with comparable quartic cou- 
plings. Right graph shows the slow-roll parameter, e, as a function of number of e-folds. 

The trajectory in <f)~x plane is given in left graph of Figj3j It contains both the analytical 
curve Eq. (jVT7j) and the exact numerical results which indicate that the slow-roll approxi- 
mation is very well valid until the end of inflation. We have also graphed the first slow-roll 
parameter 




(V.ll) 



in terms of number of e-folds in the right graph. As one can see, e changes smoothly up 
to the end of inflation. The value of quartic couplings are chosen such that the amplitude 
of density perturbations for the mode that exit 60 e-folds before the end of inflation, k^, 
is COBE normalized, i.e. Pn{km) — 2.4 x 10~ 9 . The scalar spectral index at such scales 
is n-ji(kQ ) ~ 0.947, which is slightly smaller than single A0 4 model, but still within the 2a 
level of WMAP results. 

The amplitude of correlated entropy mode is 4.4 x 10~ 17 at the current horizon scale and 
has a spectral index of ~ 0.9837. The correlation between curvature and entropy perturba- 
tions, C, is 0.228 at such scales which shows that curvature perturbations at these scales are 
partially generated through the transformation of entropy perturbations to curvature ones. 
The amplitude of tensor perturbation is Pr(^6o) — 6.33 x 10 -10 , corresponding to r ~ 0.26 
which is on the verge of being ruled out by the future experiments. The tensor spectral 
index, n?, is ~ —0.036. One should note that the consistency relation between tensor and 
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FIG. 4: Left graph shows the trajectory for two quartic potentials with hierarchy in quartic cou- 
plings. There is a transient non-slow-roll phase in the evolution of the fields. The solid and dashed 
curved, respectively, denote the analytic and numerical results. Right graph shows the slow-roll 
parameter, e, as a function of number of e-folds. 

scalar spectra for single field inflation, r = — 8ut, changes to r = —8^(1 — C 2 ) in such 



two-field models 



23], which is confirmed by our numerics. In this sense the two-field A0 4 



theory is subject to a lower tensor-to-scalar ratio r than the single field case. Hence it may 
fall in the region of n s — r plane which is allowed by WMAP5. (c.f. 24| for the violation of 
consistency relation in the context of trans-Planckian physics). 

In the other example that we will consider, there will be a transient period of fast- 
roll evolution. To realize such a scenario there should be a hierarchy between and A x . 
Inflation will not stop, nonetheless e approaches close to unity for few e-folds. Before this 
transient period, inflation is basically driven by one of the fields and after that by the other. 
In \<p A model, the spectral index depends on the number of e-folds, N e , via the relation, 
n-ji = 1 — 3/N e . Thus if the first phase lasts less than 50 e-folds, the scalar spectral index 
easily falls outside the 2a limit of WMAP5 central value for the spectral index. Therefore 
double inflation models like 18], consistent with predictions in the CMB regime, are difficult 
to be produced with two quartic potentials. Nonetheless, we investigate this possibility that 
such a non-slow-roll phase occurs toward the end of inflation. In the following example we 
have tuned the parameters such that the non-slow-roll phase occurs in the last 10 e-folds of 
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FIG. 5: The graphs show the trajectory for two symmetry-breaking cases considered in the 
manuscript. 

inflation: 

A = 2.05 x 10~ 13 , A x = 5 x 10~ 16 

0, = 21.8 M P , X i = 9-1 M P . (V.12) 

Initial values for the fields are given 70.34 e-folds before the end of inflation. The exact 
numerical trajectory is shown in the left graph of Fig. [4] by the solid curve, whereas the 
analytic trajectory, given by formula (IV. 7j) . is shown by the dashed one. The mismatch 
between the two curves is a result of slow-roll violation at the end of inflation. We have also 
included the behavior of the slow-roll parameter, e, vs. N e which shows that for few e-folds 
e increases considerably. Tuning the values of quartic couplings such that the amplitude 
of density perturbations matches the COBE normalization, one obtains n^k^) ~ 0.941. 
The amplitude of correlated entropy perturbations are 1.79 x 10~ 20 whose spectral index 
is ~ 0.958. The relative cross correlation between curvature and entropy perturbations is 
0.411 which has a blue spectral index ~ 8.4 x 10 -3 . The amplitude of gravity waves at 
Hubble scale is ~ 6.42 x 10~ 10 , i.e. r ~ 0.328, which shows that this model is completely 
excluded by the observations. 
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B. Symmetry Breaking 



The next example we study corresponds to the symmetry breaking potential where 

y = ^0 2 (0-^) 2 + ^ X 2 (0-^) 2 - (V.13) 

Using the equations of motions for and x combined with Friedmann equation one can 
show that the curved inflationary trajectory in — x fi e ld space is 

2 2 

0(0 - H) + X(X ~ the) - y ln (20 - H) - y ln(2x - // x ) = 8M 2 iV e . (V. 14) 
We consider the following values for the couplings: 



A = 2 x 1(T 15 , \i$ = 196.168M P , /j x = 36M P , A x = A^ 




(pi = 209.439 M P , Xi = 26.678 M P (V.15) 

The trajectory in the — x plane is given in the left graph of Fig. O . The theoretical 
curve Eq. ( TV. 141) and the curve obtained from the full numerical analysis coincide with each 
other which indicates that the slow-roll approximation is very well valid up to the end of 
inflation. For this model, n-ji ~ 0.963, which is compatible with WMAP 5 years results. The 
amplitude and spectral index of correlated entropy mode are respectively 1.4 x 10~ 13 and 
0.978 at horizon scales. The amplitude of correlation factor between curvature and entropy 
mode is C ~ 0.602. For gravity waves, the amplitude and spectral index are respectively, 
P T (£; 6 o) = 2.618 x 10~ 10 , i.e. r ~ 0.107, and n T ~ -0.041. 

Another interesting case that has a nice geometric interpretation is when N\ = N 2 and 
<pi > Xi- 111 the string theory picture, this corresponds to two stacks of D3-branes where in 
the background of RR six form, CVg), two of their perpendicular dimensions are blown- up to 
two concentric spheres, one of which has a radius smaller than \x and the other one bigger 
than [l. The bigger one shrinks and the smaller one expands. These two spheres collide 
when their radii reaches \x and pass through each other. This incident happens well after 
the end of inflation and during preheating and may have interesting consequences. One set 
of parameters for which, the above scenario is realized is given below: 

A^ = A x = 6 x 10~ 14 , = fi x = 26M P , 

(P t = A2M P , x* = 18M P , (V.16) 
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where the initial conditions for the fields are given approximately 73.28 e-folds before the end 
of inflation. The trajectory in the <f)—x plane is given on the right graph of Fig. [5j Fixing the 
amplitude of curvature perturbations for the modes that exit the horizon 60 e-folds before 
the end of inflation, one obtains the scalar spectral index, n-jz ~ 0.954. The amplitude 
and spectral index of correlated entropy perturbations are respectively, 4.61 x 10~ 14 and 
ns — 0.963. The relative cross correlation between the corresponding entropy mode and 
curvature perturbations is C ~ 0.78. The amplitude of tensor perturbations is 3.16 x 10~ 10 , 
i.e. r ~ 0.13. 



VI. DISCUSSIONS 

In this work we continued the analysis of the Matrix Inflation (M-flation) model we 
proposed in [4|. Our analysis was focused along two directions. We first tried to map 
the landscape of the inflationary models arising from M-flation. As discussed, due to rela- 
tive simplicity of our model compared to the theories obtained from generic string theory 
(Calabi-Yau) compactifications, we can give a complete map of the landscape, at least classi- 
cally. This enabled us to take first steps toward studying quantum effects on the landscape. 
These are basically tunneling (instanton effects in the scalar and/or gravity theories). As 
we discussed, however, such transition amplitudes are too small to bring us out of the lo- 
cal minimum. In this sense M-flation provides a toy model with a fairly rich and at the 
same time tractable landscape. This is more remarkable noting that M-flation naturally 
and quite generically occurs in string theory setting as low energy dynamics of D-branes 
in certain background fluxes. However, in calculating the potentials from the dynamics of 
coincident branes in an appropriate flux compactification, we limited ourselves to low energy 



dynamics (lowest order in a') and did not take into account the 



Dack-reactions of the com- 



pactifications and the moduli stabilization on the potential 
can have significant effects. 
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331 ] . These back-reactions 



As we discussed if we start with a field configuration for which $j are of the form (1H.5|) 



our theory (at classical level) effectively reduces to an n-field inflationary model. If at t = 
we start with the sum of two such solutions, which one leads to an n\ field model and the 
other to an n 2 field model with n\ + n 2 = N, then the dynamics of the theory is such that 
we will not simply get an n\ + n 2 field model and all 3N 2 fields will eventually be turned on. 
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It may, however, happen that at different points in time one can approximate the theory 
with an effective multi-field model, in which the number of fields may change in time. This 



provides us with a situation similar to the one discussed in 26]. 
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tjs, interesting to further 
within our M-flation 



study this line and examine the idea of "meandering inflation" 
setting which provides us with a tractable landscape. 

Next, as a show case, we considered the two-field inflationary scenarios which arise from 
M-flation and studied some of its specific features, including the entropy modes, the spectral 
index n s and the tensor-to-scalar ratio r. 

Here we did not study the preheating scenarios which naturally arise in a generic n- 
giant vacuum of M-flation. We expect the analysis of the preheating and resonant particle 
creations by non-adiabatic fluctuations of the scalar fields to be similar to the single-giant 
case studied in |4|. The n-giant case, however, is expected to have its own novel features 
too. For example, in the corresponding spherical D-brane picture, one expects an excessive 
particle production when giants (spherical branes) pass through each other at the end of 



inflation 
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35[. When the branes become coincident some of the modes of open strings 
stretched between the giants reach their minimum mass causing a resonant production of 
these modes, leading to a very effective preheating scenario. This point deserves further 
analysis which will be carried out and presented elsewhere. 
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Appendix A: Two Field Perturbation Theory 

Here we present the perturbation theory of two field models in some details. In mod- 
els with multiple inflaton fields, the field perturbations are decomposed into perturbation 
tangential to the background inflationary trajectory, the adiabatic perturbation, and the 
perturbations orthogonal to the background trajectory, the entropy perturbations. For an 
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extensive review see and the references therein. 

Following 38|, the velocity in the field space is a = \/ ip 2 + x 2 and we can define the 
polar angle in the field space as 

cos9 = Lp/a. (A.l) 
It is now useful to define the following Mukhanov- Sasaki variables: 

Q a = cos 6Q ip + sin6 Q x , (A.2) 

where 

Q v = 6<p+£$ and Q x = 5 X + £$. (A.3) 

We work in the longitudinal gauge where in the absence of any anisotropic stress-energy 
tensor the perturbed metric takes the following form [39]: 

ds 2 = -(1 + 2$(t,x))dt 2 + a(t) 2 (l - 2$(t,x))dx 2 . (A.4) 

In the flat gauge, Q a represents the field perturbations along the velocity in the field 
space. Q a is also related to the commonly used curvature perturbation, 1Z, of the comoving 
hypersurface via 

K = -Q a . (A.5) 
a 

Similarly the isocurvature perturbation is: 

5s = — sin 6 Q v + cos 6 Q x . (A. 6) 

It describes field perturbation perpendicular to the field velocity in the field space and, by 
analogy with 71, we can define a rescaled entropy perturbation, S, through 

S = —5s . (A.7) 
a 

The transformations described above basically amount to introducing a new orthonormal 
basis in the field space, defined by vectors 

E a = (E%,E%) = (cos#,sin6>) , (A.8) 
E s = (Ef, E*) = (- sin 9, cos 9) , (A.9) 

which turn out to be useful to express various derivatives of the potential with respect to 
the curvature and isocurvature perturbations. Employing an implicit summation over the 
indices I, J £ {<p, x}, one thus finds 

V a = ElVj, V S = E I S V I) (A. 10) 
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and 

V aa = ElEiV u , V as = ElEfVu , V ss = E^Vu . (A. 1 1) 

for the first and second derivatives. 

By combining the Klein-Gordon equations for the background scalar fields one obtains 
the background EOMs along the curvature and isocurvature directions 

a + 3Ha + V a = 0, (A. 12) 

9 = -—. (A.13) 
a 

With help of these equations, one can show that the EOMs for curvature and isocurvature 
perturbations become 

(k 2 \ 2V ■ 

Q a + 3HQ a + f — + C aa J Q a + -^Ss + C as Ss = 0, (A. 14) 

(k 2 \ 2V ■ 

5s + 3H5s + ( — + C ss \ 5 s - —^Q a + C sa Q a = 0, (A. 15) 

with coefficients given by 

c °° = v °°-{t) + m 2 h + W 2 -Kh 2 (A - 16) 

C as -6H- + — + 2V as + (A.17) 

C ss = V ss -(^pj (A. 18) 

C sa - ^11-- — + ^. (A.19) 



The power spectra of curvature (adiabatic), entropic perturbations and correlation spec- 
trum are defined, respectively, as 

= ^ <^*> " k ') » = ^ (52$,) 5 3 (k - k') . (A.20) 



CWsW = — (7^<S k ,> 5 3 (k - k') (A.21) 



The correlation spectrum is defined as : 

2^ 

The correlation is also often quantified in terms of relative correlation coefficient, C, which 
takes values between and 1 and indicates to what extent final curvature perturbations 
result from interactions with entropy perturbations. It is defines as 

C(k) = - lClzs{k)l (A.22) 
VVn(k)V s (k) 
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The curvature and isocurvature perturbations are evolved by assuming initially, at con- 
formal time Ti, a Bunch-Davies vacuum. Therefore, when the wavelength of the two types 
of perturbations is initially much smaller than the Hubble radius, k 3> aH, we impose the 
initial conditions 



-iki 



-ikTi 



2k 



and 



5s(Ti 



2k 



(A.23) 



Inside the horizon these two modes are independent, because their corresponding EOMs, 
eqs. (IA.14|) and fl A. 15[) . are independent in the limit k ^> aH. However, this does not hold 



when the modes leave the horizon 
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Appendix B: Phenomenologically Viable Models in Case III 

In case III, there is a nontrivial metastable giant vacuum at 0o which has an energy larger 
than the true vacuum at <p = 0. If the primordial inflation occurs in the region <p > <f) Q , 
it is possible that the inflaton gets trapped at the giant vacuum in O an d a secondary 
phase of old inflation occurs. This resurrects the graceful exit problem. However one may 
avoid this problem, if the universe can tunnel from the metastable vacuum. This could 



occur in two ways: through Hawking-Moss 



15] or Coleman-De Luccia [14j tunneling. In 



the former, the universe can exit from the false vacuum expansion through a homogeneous 
bubble solution whose radius is greater than that of de Sitter space. The phase transition 
happens simultaneously everywhere and the universe jumps from its metastable vacuum to 
the maximum of the potential and subsequently rolls downhill to its global minimum due 
to its tachyonic perturbative mode. This occurs if m 2 < 2H 2 , where m 2 = d 2 V / d(f) 2 \(f)Q. 
In the latter, which occurs when m 2 > 2H 2 , there is a solution that interpolates between 
two vacua directly. Below we will show that transition from the false to the true vacuum is 
impossible, unless the height of bump is very small. In the Hawking-Moss case, the inflaton 
rolls downhill after the transition. Taking into account the COBE normalization for the 
amplitude of density perturbations for the modes that exit the horizon during the subsequent 
slow-roll phase in the region <p < <fto, we show that Hawking-Moss phase transition is not 
possible for any value of parameters. 

First let us calculate the nucleation rate with Coleman-De Luccia phase transition. The 
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nucleation rate via the Colemand-De Luccia instanton, Tcd, is given as 

T CD = M 4 exp(-S E ) J (B.l) 

where M 2 ~ 0(V") at the false vacuum and we approximate it with V" hereafter. Se is 
the Euclidean action for the bounce solution that interpolates between the false and true 
vacua. Phase transition can occur if 2 

IF ~ ^ (B - 2) 

To estimate the Euclidean action, we approximate the potential by a triangle and will use 
the results of 2^|. Following them, we denote the false and true vacua as 0+ and 0_ with 
amplitudes V+ and V_, respectively. We also designate the local maximum of the potential 
at 4>t with amplitude Vr- If 

^V /2 > 2 *»- , (B.3) 
AV+J -A^-A^' y ' 

we then have the following expression for Euclidean action: 

3 (VTTc- 1) V AV+ / 

In the above A0± = ±(0t — 0), AV± = Vr — V± and c is the ratio of gradients of potential 
on either side of the local maximum: 

c=y> ( B - 5 ) 



where 



AV+ 

A± = -r-±. (B.6) 
A0 ± 



If condition (1B.3j) is not satisfied, then the Euclidean action is given by the following ex- 
pression 

S» = 

96 

where, 



S E = ^\ 2 + RU-Pl + 3cf3 2 + {3_ + 3c/3 + /?! - c 2 (3 3 _), (B.7) 



2 In this appendix we measure dimensionful parameters like k and m in units of Planck mass and for the 
ease of notation set Mp = 1. 
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and 

Equipped with the following formulae, we calculate the Coleman-De Luccia nucleation rate 
in the case above. Since, 

T = cos 0) (B.10) 

A^ 

+ = ^(1 + COS 0) (B.ll) 
A^ 

0_ = (B.12) 



Inequality ( TB.3j) is satisfied only when does not fall in the interval 
(arccos(|), 1.4014969565]. Inside this interval, the Euclidean action is given by ex- 
pression ( IB. 71) which yields: 

S E = £M (B.13) 

A<£ 

where /i(0) is a complicated function, but it is enough to know that it is monotonically 
decreasing as a function of and starts from infinity at arccos(|) and decreases to 344.477 
at 1.4014969565. Since A^ < 1, the nucleation rate remains negligible in this interval of 0. 
Outside this interval, the Euclidean action is given by expression flB.4j) : 



S E = (B.14) 



where 

1287r 2 (l + 3 cos3 Q~ 5 cos2 Q+cos e+i 

4 • 



/ 2 (e) = 8cos2Q — J — 4 . (B.i5) 



i I 3 cos 3 Q— 5 cos 2 Q+cos €>+! \ -r 

1 + 8cos 2 e I l 



/ 2 (0) reaches zero when tends to |. Close to |, /2(@) behaves like 

/ 2 (0) 1024tt 2 



0-^ 3 
2. 



(B.16) 



Mass parameter and Hubble parameter at + can be calculated and they turn out to be as 
follows: 



2k cos 0(1 + cos 0) 

A^ 2 ~m(0+) 2 = J (B.17) 

A^ 

g2(0+) = ^(l + cosem-3cos9) (Blg) 
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Plugging these expressions into eq. ( IB.2j) and expanding its L.H.S. around = |, one 
obtains 

Too 5184A{(9 - |) 2 ^ 1024 T 2 (e-f)^ 

a exp ( s; J ■ (B19) 



If we call x = — 102477 ( e ^Z 2 ) ; eq.1 1B.2l) takes the form 



81A 14/3 Q 
_J_ x Vs eM _ x) >_ (B , 0) 

Noting that < Q < n/2, x varies between and 12 ^ 7r . This equation has solutions, if the 
maximum of L.H.S. of equation is bigger or equal to 9/4it. Since the maximum of L.H.S. of 
equation happens at x — 2/3, we obtain the following constraints on the parameters: 

-V > exp(2/3) ~ 0.1648. (B.21) 

k\ 36 

If this condition is satisfied, since the L.H.S. of the equation has a Maxwellian shape, in 
general we will have two solutions X\ and X2 for the equation which correspond to angles 0i 
and 6 2 where § - 4 = (i^) V3 . For 

^V /3 <--e<P^V /3 , (B.22) 
1024vr 2 J ~ 2 - Vl024vr 2 / ' K J 

the Coleman-De Luccia nucleation rate is bigger than the critical rate and inflation can end 
via Coleman-De Luccia first order phase transition. On the other hand, Coleman-De Luccia 
solution only exist when m 2 > 2H 2 . This results in addition constraint on the couplings 

4 36 cos m 2 on 

A 2 (1 + cos0) 2 (1-3cos0)' 1 ' ' 



or since is very close to 7r/2 



7T „ f\ 



< — ^ (B.24) 
2 36A 2 V ; 

If the parameters and A^, are tuned such that satisfies (IB. 241) and (1B.22j) simultaneously, 

we will have successful Coleman De-Luccia phase transition at the end of slow-roll inflation. 

The potential expanded around = </>+, with ~ |, is: 

"to = + + ^. (B-25) 

where if = <p — <p + . For ip's bigger than but close to zero 77, the second slow- roll parameter, 
is bigger than one. However for (p ~ O(10)Mp, the potential can sustain chaotic inflation. 

29 



It is possible to show that one can satisfy the above constraints and the ones from WMAP 
simultaneously. For example consider the following values for and k$ 

A = 1.7 x 1(T 18 , = 2.5 x 1(T 21 M P 

The constraint (IB.21I) is clearly satisfied and one finds the following allowed interval from 
combining flR22l) and (1R24D : 

2.694 x 10~ 8 < | - 6 < 6.007 x 10~ 8 . (B.26) 

For these values of couplings the amplitude of density perturbations 60 e-folds before the 
end of inflation matches the COBE normalization, 2 x 10 -5 and the spectral index is 0.9508 
which is within the WMAP5 2a level. 

We also would like to show that the phase transition from the false vacuum to the true one 
cannot happen via Hawking- Moss phase transition, if one demands to produce the observed 
amplitude of density perturbations in the subsequent slow-roll phase. The nucleation rate 
for Hawking-Moss phase transition is given by an expression similar to ( IB.lj) where now the 
Euclidean action is calculated for the Euclidean solution that interpolates between <b + and 

n 

4>t- According to [15(, the Euclidean action is given by 

3eu = (B.27) 

or 

1536tt 2 AJ cos 3 6 

Shm = 4 • 6rvi o nvi , o rvr (B.28) 

k^l sin 6(1 — 3 cos 6) (1 + 3 cos 6) 

Shm approaches zero as we approach 6 = |. This justifies that we expand it around this 

point: 



1536tt 2 A 



2\3 



7T 



Using (11117]) & ( jBTgj) , we have 



Shm = ^ 6-- +■■• (B.29) 



TT\2 



; , 5184A4(6-f ; 
Thm/H* w ^ ^exp 



1536vr 2 Ai / n\Z 



Kl ■ /cl V 2 



(B.30) 



We now follow the same procedure we did to analyze the behavior of the nucleation rate in 
the case of Coleman-De Lucica phase transition. We define the variable x in the following 
manner 



15367r 2 A 3 / 7T\3 

X " ^( -2^ ' (R31) 
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In terms of which the equation Thm/H 4 > 9/4ir takes the form 

3 1C V 3 A 2 q 

jj= * x 2/3 exp(-x) > f-. (B.32) 

In order to have a non-empty set of solution for this equation, the couplings have to satisfy 
the condition 

A3, 9 / 7T 2 \ 1/3 

This condition on the parameters is necessary to have successful phase Hawking- Moss phase 
transition from <p + to <px- If one requires that the inflationary scenario that arises in the 
region <p < <f)j- be phenomenologically viable too, one should also take into account the 
constraint that amplitude of density perturbations sets on the parameters (1B.21j) . The 
amplitude of density perturbations 60 e-folds before the end of inflation is set by the ratio 

m ' 

S H ~ —^±N 2 e ~ 2 x 10" 5 (B.34) 
5tt m 

where N e is set equal to 60 which corresponds to the moment CMB scales have left the 
horizon. This equation sets the following relation between m and A 

m ~ 2.29 x 10%, (B.35) 

Since G is very close to |, the condition necessary to have inflection point inflation, = 
m \f^>i s t m holds to a very good approximation, i.e. 

~ 2.29 x 10 7 A 3/2 . (B.36) 



Plugging this condition into inequality (IB. 331) . one obtains the L.H.S. to be just a number, i.e. 
1.53 x 10~ 10 , which is much smaller than the R.H.S. of inequality. This means that Hawking- 
Moss phase transition with a subsequent inflationary period with the correct amplitude of 



density perturbations is not possible. This might be problematic for MSSM inflation [421 ]. 
since it shows that in the case one deviates slightly from inflection point condition such 
that the potential acquires a metastable minimum at 0o 7^ 0, the probability of Hawking- 
Moss tunneling from the false vacuum to the maximum of the potential is small, if one 
requires to obtain the correct amplitude for density perturbations in the subsequent slow- 



roll phase. This is contrary to the result of 43j, which had claimed otherwise. In fact, 43] 



had only proven the existence of Euclidean instantonic solution that extrapola tes between 
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the metastable minimum and maximum of the potential and had not calculated the rate by 
which this tunneling occurs. Our computation, comparing the nucleation and the expansion 
rates, excludes such a possibility. 
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